Radiative coupling and weak lasing of exciton-polariton condensates 
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In spite of having finite life-time exciton-polaritons in microcavities are known to condense at 
strong enough pumping of the reservoir. We present an analytical theory of such Bose-condensates 
on a set of localized one-particle states: condensation centers. To understand physics of these arrays 
one has to supplement the Josephson coupling by the radiative coupling caused by the interference of 
the light emitted by different centers. Combination of these couplings with the one-site interaction 
between the bosons leads to a rich nonlinear dynamics. In particular, a new regime of radiation 
appears. We call it weak lasing: the centers have macroscopic occupations and radiate coherently, 
but the coupling alone is sufficient for stabilization. The system can have several stable states and 
switch between them. Moreover, the time reversal symmetry in this regime is, as a rule, broken. 

PACS numbers: 67.10.-j,42.25.Kb, 78.67.-n, 71.36. +c 



Introduction. — Condensation of exciton-polaritons 
(EP) in semiconductor microcavities formed by two 
Bragg mirrors with a quantum well between them was 
recently discovered [lH5|- As experimental implemen- 
tation of the Bose-Einstein condensation (BEC) these 
systems have some advantages as compared with cold 
gases, e.g., the vortex dynamics and the superfluid 
motion can be accessed by optical methods [6|-l8|. At 
the same time, in contrast with the atomic BEC the 
EP-condensates are not in the thermodynamic equilib- 
rium. To reach a macroscopic number of EPs, which 
life-time is finite, one needs an outside pumping. Being 
driven, the EP-condensates differ fundamentally from 
the conventional BEC-systems. In particular, EPs can 
condense into one-particle excited state d, Q or even 
into several excited states 13, 



It is not unusual for a one-particle bosonic state to be 
localized. Such states either formed by disorder or inten- 
tionally prepared 0, Q can serve as condensation centers 
(CC). The bosons (photons, excitons, polaritons, etc.) 
arrive to each CC from an incoherent reservoir created 
by the pumping and escape in the form of light radiation. 
At low enough pumping one should expect a system of 
disconnected BEC droplets emitting light of different fre- 
quencies. As the pumping increases these sources of ra- 
diation tend to become coherent and to shine as a laser. 
This effect resembles the Insulator-Superfluid transition, 
but the similarity is limited by the dissipative nature of 
the system. In this paper we develop the theory of lasing 
in non-equilibrium BEC, which allows one to understand 
qualitatively and describe quantitatively a number of un- 
usual and surprising experimental observations. 

Existing theoretical approach to non-equilibrium EP 
condensates [l^ . [isj is based on the Gross-Pitaevskii 
equation (GPE) modified to account for the finite life- 
time of the bosons, the continuous feed of the condensate 
from the reservoir, and effects of reservoir depletion. Nu- 
merical analysis seems to agree qualitatively with some 
experiments. However, the nonlinear GPE can have sev- 



eral distinct solutions and the choice is ambiguous. 

Our starting idea is that the sensitivity of the bosonic 
life-time to the symmetry of the wave- function dominates 
this choice. If two CCs are close to the resonance and not 
too far apart, the interference of the light emitted from 
different CCs is constructive for symmetric (bonding) 
state and destructive for antisymmetric (anti-bonding) 
state, i.e., the bosons live longer in the anti-bonding case. 
Existence of several states with different life-times leads 
to a crossover range of pumping strengths (rather than 
a single threshold) where the income and outcome rates 
match. Within this range a coherent condensed state is 
formed. In this respect, the non-equilibrium EP conden- 
sation in disordered cavities resembles the random lasing 
phenomena: the lasing states also possess the longest 
life-times. Important difference is the nonlinearity: in- 
teraction between the EPs is able to synchronize the fre- 
quencies of CCs [ij, ll5|, which have no reason to be in 
resonance at low occupancies. One can call this regime 
"weak lasing": the occupations of the CCs are macro- 
scopic and the radiation is coherent, but the state is sta- 
bilized by the coupling between CCs rather than, e.g., by 
the depletion of the reservoir. 

Below we show that CCs indeed become coherent and 
form a particular long-living condensate. In certain range 
of the parameters the system can have several linearly 
stable states, and switch between them. We argue that 
the properties of these states and the switches naturally 
explain otherwise mysterious experimental observations. 

Formalism. — In the familiar Glauber-Sudarshan for- 
malism of coherent states 16] the occupation n and 
the phase (p of each CC are encoded in the complex 
number z — ^/ne"^. N CCs arc thus characterized 
by A^-component complex vector Z = {zi, Z2, . . . , zjv}, 
where corresponds to the CC number /i. Each state 
of this system is fully described by the density matrix 
p(Z, Z* ,t), which evolves in time according to (ft = 1) 
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Here is the incoming rate of the bosons to the /x-th 
CC, = d/dZfi, dot indicates the time derivative, and 
% is the complex Hamiltonian function 

The Hamiltonian function of an isolated CC is 

H(i)(|z^n = i (r. - W,) \z,\^ + iH,{\z,\% (3) 

where is the escape rate, H^^{n^i) is the energy of 
bosons in the /i-th CC, which alone emits light with the 
frequency = dH^^(jif^)/dn^j_. For weak interaction 



+ \a.f^n^_,, (4) 



with one-particle energy and the coupling constant 
> 0. The coupling is weak, ^ w^, and £7^ slightly 
increases with the pumping. In general, the income rates 
are dynamical variables determined by the occupa- 
tions of the reservoir states. However, in (1-3) enter 
as Z-independent external parameters. We adopt this 
simplifying assumption, which is valid when the life-time 
of the reservoir excitations is not limited by the inelas- 
tic relaxation into CCs. It is quite straightforward to 
generalize our approach to include the dynamics of . 

The coupling V^j, between CCs in ^ is bilinear, i.e., 
the bosons from different CCs do not interact, and it 
consists of radiative (dissipative) and Josephson (non- 
dissipative) parts, V^j^ = 7^1/ — iJ^jiy- The Josephson cou- 
pling favors symmetric ground state, i.e., J ^ 0. The 
radiative coupling is due to the interference of the 
light from different CCs. The matrices 7^1^ and J^^ are 
Hermitian and have zero diagonal matrix elements, which 
are absorbed by the first term in Besides, we assume 
no explicit time reversal symmetry violation and thus 
7^1/ and Jfiu are real. Note that the Josephson coupling 
caused by the tunneling is exponentially small, when CCs 
are spatially separated, while 7^^ can be still substantial. 

The time evolution of the density matrix without 
Josephson coupling should reflect the following physi- 
cal picture. According to ([2]) and ([3]) both 3?z^ and 
QZfi oscillate with the high frequency £7^ . As long as 
flfj, ^ Slj^, coupling between the CCs /i and ly averages 
out as time exceeds jJl^ — ri,y|~^. However the interac- 
tion between the bosons synchronizes the frequencies. As 
soon as this happens, the interference suppresses the es- 
cape if the phases of CCs are opposite. This suppression 
in turn stabilizes this nontrivial stationary state (NTSS) 
with high occupations and coherent radiation of the CCs. 
NTSS differs drastically from the trivial stationary state 
(TSS) with small and different emission frequencies. 
We emphasize that NTSSs appear due to the off-diagonal 
radiative coupling even at J^^, = 0. 

At this point we adopt the Langevin approach, which 
remains valid for a complex Hamiltonian function. We 



start with the formal solution of Eq. ([Ij for p(Z,Z*,i) 
with the initial condition poiZ, Z*) in the form of Feyn- 
man path integral over the "trajectories" Z{t) 

p(Z,Z*,i)= /po(Zo,Z5)d2^Zor'*^ ^I?Z(T)e-^ (5) 

J JZ(0)=Zo 
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Note that the pumping intensity of the /i-th CC plays 
the role of its effective temperature. 

One can rewrite ([5]) as the average ((5(Z — Z*^-''-'(t))), 

where z\I\t) are the solutions of equations 
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with ffj,{t) being a realization of the Gaussian random 
processes with zero mean {fp,{t)) = and (5-like two-point 
correlation function {ff_i{t)f*{t')) = W^6^i,S{t - t'). 

Our analysis of ([7]) consists of the following steps. First 
we omit noise to find the stationary points. Next we 
analyze the stability of these points and thus identify the 
long-living stationary states. If such states coexist, the 
time the system spends in each basin of attraction is finite 
due to the noise /^(t). The set of these times determines 
the behavior of the system. Below we evaluate these 
times for the simplest case. The final step — the analysis 
of the fluctuations of Z in the vicinity of the stationary 
points, which determine the line-shapes of the radiation 
and its statistics — lies beyond the scope of this letter and 
will be discussed elsewhere. 

Substitution of (2,3) into Q at = leads to 
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Being similar to the equations used in 1^ Eqs. ([5]) dif- 
fer in two aspects. First one is the radiative coupling. 
The second difference is the absence of the reservoir dy- 
namics. We believe that the depletion of reservoir is not 
qualitatively important in the weak lasing regime. 

Two Condensation Centers. — From now on we restrict 
ourselves by the case of two coupled CCs. The density 
matrix p^'^^{zi,Z2) depends on four real variables, how- 
ever, the total phase, is irrelevant to what follows. We 
choose to parameterize the remaining degrees of free- 
dom by three components of a 3D classical spin S = 
{Sx, Sy,Sz) using the Pauli matrices cr^y. 
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= (|zip-|-|z2p)^/4. For simplicity 
Wi) = (F2 - W2) = 5; a small 
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difference between (Fi — Wi) and — W2) would not 
affect the results qualitatively. From ([8]) and ([9]) we have 



-SV{^H) - [S X V(Q?-H)] 



(10) 



where V = d/dS. Calling uj ^ oji — uj2, a = {ai +a2)/2, 
and (3 = (02 — cti)/2, we rewrite (fTOl) explicitly 



[oj + aS,- (3S]Sy, (11a) 
-[oj + aS,-f3S]S,, (lib) 
(11c) 

- jS^. We used 7^^ = jaf^^ 



Sx = -gSx - -fS ~ 

Sz — gSz J Sy, 

and accordingly S = —gS ^ 
and Jf_,i, ^ Ja^^. 

The stationary states are solutions of (|lll) at S = 0. 
It turns out that, apart from TSS S = 0, as long as < 
5 < 7 there exist NTSSs. In polar coordinates (S, ip) 
their positions are 
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— arctan r, 



cosi9 = —Jr/j, 
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Since r = ±|r|, Eq. (IT^) for S gives two solutions. 
Each positive solution corresponds to a NTSS. S changes 
the sign together with either its numerator or its denom- 
inator, i.e., g has two bifurcation values, go and 500: 



^'ffo = (^'+.9o)(7'-3o), (13a) 
{f3^j' + a'j')gl^j'j\a^~n (13b) 



At g = go one of the NTSSs merges with the TSS. At 
g = ffoo and r < the value of S diverges. 

To manifest itself as an attractor the NTSS has to be 
linearly stable. The stability diagram Fig. 1 follows from 
the analysis of the Lyapunov exponents. The TSS loses 
its stability at gg. The stability of NTSS with a given S 
()12b ) changes at g — gc defined by the equation 

[( J2 - g^)l3jr + (72 + g^)aj] gS = 27(5^ + 7^^). 

(14) 

Here we are not aiming to describe the nonlinear dy- 
namics of the two CCs system in all details. In particular, 
we discussed only the stationary states. However, there 
are numerical evidences for the existence of the limiting 
cycle (LC) in the regions IV and V in Fig. 1. Moreover, 
within the region V all of the stationary states are unsta- 
ble, i.e., there should be at least one stable LC. Analysis 
of stable time- dependent solutions and their experimental 
manifestations is a subject for further studies. 

Kramers transitions. — In the regions II and VII in 
Fig.l two stable stationary states coexist, which leads 
to an additional problem: the noise term in ^ causes 
switching between these states. The properties of the ra- 
diation from the two states are very different. If both 
escape rates are small enough, this switching can be 
observed directly. Otherwise we expect a combination 
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FIG. 1. Bifurcation and stability curves, according to Eqs. 
(O and (O, for J/7 = 0.25 and a = 2/3. Stable(unstable) 
TSS and NTSS are indicated by STR(UTR) and SNT(UNT) 
in each domain. Zoomed domain V is shown in the insert. 



of the two signals, the weights being determined by the 
relation between the escape rates. 

One can evaluate r^^ analytically when J — and 
according to ([13b ) 500 = 0. In this case the region HI 
disappears and the bistable region II is defined by the 
inequalities uj > -f > g > 0. At g = the system is 
conservative and integrable, i.e., pT|) has two integrals 
of motion: [see (HB)] and A = 2P{-fSy+ijS) - 13^3"^. 
Stationary states correspond to 6*2 = for all g > 0. 
For TSS A = 0, while for NTSSs A = (w ± 7)^ = A± 
(the stable and unstable NTSS correspond to -I- and — , 
respectively). Other values oi A a.t g — label periodic 
trajectories (see Fig. 2). Slow time evolution of A at 
small g > is governed by 
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+ F{t), {F{t)F{t'))=f}TS{t-t'), (15) 
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dE 
du) 



(16) 



where E = K+i - 27ruj9{A^ - A), 
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R^^/^0{R)d(p, R= {lu+ J sin ip)^- A, (17) 



with e{x > 0) = 1 and 9{x < 0) = 0. 

Equation (|15p is a typical case of the one-dimensional 
Kramers problem — classical transitions between two po- 
tential minima due to the thermal noise — slightly modi- 
fied by the fact that "effective temperature" pB)) depends 
on the "coordinate" A. Using the familiar solution 

[13- 

[19I I we estimate the escape rates from the stable TSS and 
NTSS as 
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FIG. 2. Flow near NTSSs and TSS for J = and u) = 67. 
(a) The saddle trajectories for g — O.257. (b) Several closed 
trajectories for g = 0. 



where = A_ = {uj — 7)^, Ag = and Ag = A+ = 
(a; + 7)2 for TSS and NTSS, respectively. This approach 
is vaHd only provided that the exponential factor is small, 
i.e., when w//? >• 1 and thus the occupations at NTSS 
are large. 

One can use (fTS]) to compare the steady state occupa- 
tions of TSS and NTSS. It turns out that for small cj the 
NTSS prevails. The escape rate from NTSS is approxi- 
mately independent of w, while the escape rate from TSS 
decreases with uj. The rates become equal at w ~ 4.I7. 
For bigger uj TSS dominates. 

Comparison to experimental observations. — We believe 
that a number of recently observed features of the radi- 
ation of the disordered exciton-polariton structures can 
be naturally interpreted in the framework of our theory. 
Large, 7r/2 < ip < 37r/2, phase differences between the 
CCs, which manifests itself in the dip of the radiation in- 
tensity /(k^) at small transverse wave- numbers kj^ was 
reported by several groups ^, iQi-jllj . A detailed analy- 
sis of the radiation from disordered CdTe structures was 
performed in 



Hi 



It was found that (i) the number of 
CCs differs from the number of radiated frequencies; (ii) 
some frequencies are radiated from several CCs, and (iii) 
its distribution in the k^-space as a rule has a charac- 
teristic annular shape. The last but not the least impor- 
tant observation in [ll| is (iv) the clear absence of the 
k_L -k± symmetry, I{k±) ^ I{-'k±). 

Features (i), (ii) is exactly what we should expect: for 
example, as long as TSS and NTSS in the system of two 
CCs coexist each CC radiates its own frequency and in 
addition they together radiate the third one, i.e. two 
CCs radiate three frequencies and one of them is radiated 
from both of them. The number of collectively radiated 
frequencies quickly increases with the number of CCs. 
(iii) had already been discussed - it reflects the large 
phase differences favored by the radiative coupling. 

(iv) The time inversion symmetry implies that the 
phase differences between different CCs can be either 
or TT. One can check that under this condition /(kx) — 
/(-k_L). However, according to (HDi) in the NTSS the 



phase mismatch = ^2 — 0i between the two CCs is nei- 
ther nor TT. We believe that by observing states with 
broken time inversion symmetry authors of ll| experi- 
mentally proved existence of the radiative coupling. 

In conclusion, we considered radiation of coupled cen- 
ters of condensation of exciton-polaritons in the weak las- 
ing regime, where the occupations ni and n2 of the cen- 
ters are macroscopically large, while the occupations of 
the states of the incoherent reservoir remain independent 
on rti_2. Apart from the usual Josephson coupling be- 
tween the centers we took into account the radiative cou- 
pling due to the interference of the two sources of light. 
The radiative coupling turns out to be responsible for the 
very existence of the weak lasing regime and crucial for 
understanding the experiments on coupled condensates. 
The onsite interaction between the bosons makes the dy- 
namics of the system nonlinear and rich. In particular, 
depending on the parameters there can be more than one 
stable stationary states, and some of them are nontrivial: 
in contrast with the trivial stationary state (independent 
sources of radiation) the radiation of the two centers is 
synchronized in frequency. The phase difference is locked 
thus breaking the time-reversal symmetry. This symme- 
try breaking was recently observed 
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